A limit of the confluent Heun equation and the Schrodinger equation for an inverted 

potential and for an electric dipole 
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We reexamine and extend a group of solutions in series of Bessel functions for a limiting case of 
the confluent Heun equation and, then, apply such solutions to the one- dimensional Schrodinger 
equation with an inverted quasi-exactly solvable potential as well as to the angular equation for an 
electron in the field of a point electric dipole. For the first problem we find finite- and infinite-series 
solutions which are convergent and bounded for any value of the independent variable. For the 
angular equation, we also find expansions in series of Jacobi polynomials. 



I. INTRODUCTION 

Firstly we revise and extend some solutions in series 
of Bessel functions for a limiting case of the confluent 
Heun equation (CHE) - the CHE is also called gener- 
alized spheroidal wave equation [H, [3j @1- Then we ap- 
ply these solutions to the one-dimensional Schrodinger 
equation with an inverted quasi-exactly solvable poten- 
tial considered by Cho and Ho [1, [H], as well as to the 
angular equation for an electron in the field of a point 
electric dipole 0. 

For the inverted potential we find even and odd eigen- 
states given by finite series of Bessel functions corre- 
sponding to the same energy eigenvalues, for any value of 
the parameter which characterizes the quasi-exact solv- 
ability. We obtain as well infinite-series solutions which 
are convergent and bounded over the entire range of the 
independent variable. 

For the angular equation, in addition to expansions in 
series of Bessel functions, we get solutions in series of 
Jacobi polynomials, equivalent to the ones proposed by 
Alhaidari [3]. The latter type of solutions is inapplica- 
ble to the inverted potential and this is the reason for 
regarding expansions in series of Bessel functions. 

The stationary one-dimensional Schrodinger equation 
for a particle with mass M and energy E is written as 



with 



(Pip 



£ = 



[£ - V{u)]ij = 



(1) 



2ME 



V{u) = 



2MV{x) 



where a is a constant with inverse-length dimension, h is 
the Plank constant divided by 27r, a; is the spatial coor- 
dinate and V{x) is the potential. For the case considered 
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by Cho and Ho 



(1/4) 
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V{u) = — — sinh u „ , 

4 cosh u 

[li e (— cx), cx)), l — l,2,3,---] 



(2) 



where b is a positive real constant. This is a bottomless 
potential in the sense that V(u) — > — c» when u —^ ±(X). 
If b^ < 4(/ + 1)2 - 1, it is an inverted double well; if 
b^ > 4(Z+1)2 — 1, it resembles the potential of an inverted 
oscillator. Then, we have the equation 



(ftp 



b^ 

4 



sinh M + / 



cosh^ - 



V = 0, (3) 



which, for the previous values assigned to /, is a quasi- 
exactly solvable equation. 

A quantum mechanical problem is quasi-exactly solv- 
able (QES) if one part of its energy spectrum and the 
respective eigenfunctions can be computed explicitly by 
algebraic methods [1, [^, From another viewpoint 

[ll| , a problem is QES if admits solutions given by finite 
series whose coefficients necessarily satisfy three-term or 
higher order recurrence relations, and it is exactly solv- 
able if its solutions are given by hypergeometric func- 
tions (two-term recurrence relations). This is suitable 
for problems involving Heun equations because in gen- 
eral these present finite-series solutions [l2|; they admit 
as well infinite-series solutions which could afford the rest 
of the spectrum. 

As to the second problem, the Schrodinger equation 
for an electron in the field of a point electric dipole is 
used for modeling the scattering of negative ions by polar 
molecules @. When it is separated in spherical coordi- 
nates (r, 9, if), the ^-dependence satisfies the equation 

2 1 

777 

+ C-I3COS0 ^ 6(6*) = 0, 

sin 6 

(4) 

where C is a separation constant to be determined, /3 
is the dipole moment parameter and m is the angular 
momentum conjugate to ip. 

In sections HI and IV, respectively, equations ^ 
and (j4]) are transformed into a confluent Heun equation 
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(CHE) having the form 



.(.-.o)^ + (5i+5..)f + 
[S3 + q{z - zo)] C/ - 0, (g ^ 0) 



(5) 



where zq, i?i and g are constants (if 9 = 0, the equation 
degenerates into a hypergeometric equation). This equa- 
tion is obtained by applying the so-called Whittaker-Ince 
limit [ll,[l3 

w — > 0, 1] 00 such that 2r/a; — — (6) 
to the following form of the confluent Hcun equation 3 



dU 



z{z - zo)-^ + {Bi + B2Z)— + 
[B3 ~ 27]uj{z - zo) + uj^z{z - zo)] C/ = 0, (7) 

where zq, Bt, rj and u are constants. In both cases, 
z — and z = zq are regular singular points having 
indicial exponents (0, I + Bi/zq) and (0, 1 — B2 — Bi/zo), 
respectively. However, equations ^ and ([7]) differ by 
the behavior of their solutions at the irregular singular- 
ity z = 00, namely p^ . 

U{z) ~ e±2^v/^z(i/4)-(^^/2) for Eq. ®, 



C/(z) - e^i'^-z'^iV-(B2/2) fQj. Q_ 



(8) 



We shall refer to equation ([5|) as Whittaker-Ince 

limit of the CHE ([7]). This equation also appears in 
the separation of the variables for the Laplace-Beltrami 
operator in an Eguchi-Hanson space [H, [lB| and for the 
Schrodingcr equation for an ion in the field of a elec- 
tric quadrupole in two dimensions [iTj . In addition, it 
includes the Mathieu equation as a particular case and, 
when Zo ^ 0, it leads to a double-confluent Heun equa- 
tion which arises, for example, in the scattering of ions 
by polarizable targ ets [T3 | . 

In reference the limit (O is used to flnd solutions 
for equation ([5]) from solutions of the CHE ([7]). However, 
those solutions are not sufficient to handle the previous 
problems. For this reason, from known solutions in series 
of Hankel functions for equation ([5]) , in section II we con- 
struct solutions in scries of Besscl functions of the first 
and second kind. We find that the expansions in terms of 
the functions of the first kind, under certain conditions, 
are convergent and bounded for any value of the indepen- 
dent variable. The other expansions in general converge 
for |z| > |zo| or |z — zo| > |zo|, but under special condi- 
tions converge also at | z | = | zo | or | z — zo | = | zo | . 

In section III, using the expansions in series of Besscl 
functions of the first kind, we obtain even and odd cigen- 
functions for the Cho-Ho equation ([3]), given by finite 
series and by infinite series as well. The three-term re- 
currence relations for the coefficients of the finite series 
imply that even and odd states are degenerate. A proof 



that this degeneracy takes place for any value of the pa- 
rameter I is possible because now we have a general ex- 
pression for the eingenstates. Expansions in series of the 
other Bessel functions also lead to solutions bounded for 
any value of z, but deprived of the parity property. 

In section IV we apply the solutions in series of Bessel 
functions for the angular equation of the point dipole. 
For this problem, we also get solutions in series of Jacobi 
or associated Legendrc polynomials and point out why 
this type of expansions docs not work for the inverted 
potential. 

Conclusions are reported in section V, while Appendix 
A contains some formulas used throughout the paper and 
Appendix B exhibits solutions which have been omitted 
in section II but arc used in sections HI and IV. 



II. SOLUTIONS OF THE WHITTAKER-INCE 
LIMIT OF THE CHE 

All the solutions we will use for the Whittaker-Ince 
limit of the CHE may be generated by applying the limits 
([6]) on solutions of the CHE (O, but may also be derived 
directly from ([5]). In this section we deal with a group 
of solutions given by expansions in series of four Bessel 



functions which are denoted by zi^\x) {j 
according as 



1,2,3,4), 



Zi3)(x)=i7W(x), Z(^\x) = H(^\x), 



(9) 



where Jq, and Yq, are the Bessel functions of the first and 
second kind, respectively, whereas Ha^ and Ha^ are the 
first and the second Hankel functions (a is called the 
order of the functions). Thus, a first set of solutions has 
the form 



U[^\z) = z(i-^^)/^ E(-l)"^" ^^2^"U-i (2V^) ,(10) 

[S2 ^0,-l,-2,...]. 



n=0 



where the presence of the same coefficients b'iP is due to 
the fact that the four Bessel functions satisfy the same 
difference and differential equations [l^ . These are called 
one-sided series because the summation runs from zero 
to positive infinity. 

In a region where the four expansions are valid, only 
two of them are independent. For instance, from equa- 
tions (|A.2P and (jA.SP we get 



u['\z)+ur>{z) 



(4), 



C/f \z) - C/f \z) 



(11) 
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up to a multiplicative constant. However, each expansion 
presents different behaviors at z = oo as inferred from 
equation (|A.4|) . 

In section II. A we explain how the above set of solu- 
tions has been obtained, write the general forms of the 
thrcc-tcrm recurrence relations for the series coefficients 
and present the transformations used to generate addi- 
tional sets of solutions. In section II. B we analyse some 
properties of the solutions and discuss the case of the 
Mathieu equation, while in section II. C we study the con- 
vergence of the solutions. 



A. Recurrence relations and transformations of the 
equation 

The expansions in series of Hankcl functions which ap- 
pear in (fTO|) are taken from Rcf. ^] where they have 
been expressed by series of the modified Bessel functions 

■f^2n+B2-i(±2i^/9z) as 

oo 
n=0 

[B2^0,-l,-2,---]. 
In terms of Hankel functions [1^ , 

A-2„+B.-i(-2*V9^) cx (-l)"ff<,V+i3.-i(2%/9^)' 

i^2n+B.-l(2*V^) (X (-l)"if(^VB.-l(2v^), 

where the proportionality constants do not depend on 
the summation index n. Thus, the preceding solution 
acquires the form of the solutions t/p-* and u[^\ The 
new solutions u[^^ and f/}^^ result from the properties of 
the Bessel functions, as aforementioned, and are used in 
sections III and IV. 

The recurrence relations for series coefficients 6„ and 
the respective characteristic equations (in terms of con- 
tinued fractions) assume one of the following forms [3] 

aobi + Poho = 0, | 
a„6„+i -I- /3„6„ + jnbn-i (n > 1) J 



ao7i ai72 0273 



(12) 



Pi- P2- P3- 

aoh + Poba = 0, 
ai&2 + Pibi + [a-i + 7i] bo = 0, 

Unbn+l + Pnbn + Inbn-l = (n > 2) 

ao -f 7i] ai72 "273 
= 3 z Z ' i^^i 



/3i 



P2- P3- 



ao&i + [Po + a-i] bo = 0, 

an&n+l + Pnbn + Inbn-l = (n > 1) 



where the coefficients a„, /?„ and 7„ depend on the pa- 
rameters of the equation and on the summation index n 
of the series. The second and the third types of recur- 
rence relations occur only when the parameters of the 
equation imply that a^i ^ 0, as in the case of the Math- 
ieu equation. 

If 7„ = for some n = + 1, where A^ is a natural 
number or zero, the one-sided series terminate at n = A^ 
giving finite-series solutions with < n < A^ [l^. In 
particular, for A^ = the series presents just the first 
term except when the recurrence relations are given by 
(fO)) : in this case the series presents only the first term 
if a_i -|- 7i = 0. Solutions given by finite series are also 
called quasi-polynomial solutions or Heun polynomials. 
On the other hand, the previous characteristic equations 
were established by assuming that the summation begins 
at 71 = 0. However, the series begins at n = A^ + 1 
if Un = for some n = N . In this case we must set 
n = TO + A^ + 1 and relabel the series coefficients in order 
to get series beginning at to = 0. 

The fact that in each series all the Bessel functions 
must be independent imposes certain restrictions on the 
parameters of the differential equation as in the pre- 
vious set of solutions where B2 ^ 0, — 1,— 2,--- (zln 
and ^ii^i general are proportional to each other if 
m is integer). Such restrictions also assure that the 
coefficients of the recurrence relations are well defined 
in the sense that there are no vanishing denominators. 
By transformations of variables which preserve the form 
of equation ([5]) we find new sets of solutions demand- 
ing different restrictions on the parameters. In effect, if 
U (z) = U{Bi, B2,B^]Zo, q; z) denotes one solution (or set 
of solutions) for equation ([5]), then the transformations 
=3^ and =3^, given by 

.^iU{z) = C/(Ci, C2, C3; zo, q; z), 



32U{z) - (z - zof-'''-^U{Bi,D2.D^- zo, g; z) 



^zU{z) = U{-Bi - i?2Zo, -62,-83 - qzo\zo,-q\ zq - z), 

(15) 

generate a group having eight (sets of) solutions. The 
constants Ci and Di are defined by 

2B^ 



Ci = -Bi-2zo, C2 = 2 + B2 + — , 



C3 — -B3 



D2 = 2-B2- 



Zo J \ Zo 



2Bi 
Zo 



D3 = B3 + — 
Zo 



^+B2-1 
Zo 



(16) 



Po + Oi-i = 



«o7i ai72 ^273 
Pi- P2- P3- 



(14) 



The application of these rules is straightforward but leads 
to solutions with different properties. 
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B. Two subgroups of solutions in series of Bessel 
functions 

Given the initial set of solutions, Ui , the rules and 
^2 are used as 



U['\z) ^ U^"iz) A Ul''{z) A K^'iz) (17) 

in order to generate a subgroup constituted by four sets 
of solutions. To get a second subgroup of solutions, first 
we take U^''\z) = 3^^u'f\z) and, after this, we use the 
rules 5^1 and 5^ in the same order as in (H?]), that is. 



.^1 



Ui'\z) ^ U^^'iz) A UY'iz) A U','>iz) (18) 

Then, we can check the following correspondence between 
the solutions of the two subgroups 



.^1 ttU) 



ST', 



if), 



U^\z) 



U^^^iz), (* = 1,2,3,4) 



in the sense that formally we have 



h^n r'n ' 



"n ln+1 



a. 



(»+4) (»+4) 
In+1 



(19a) 



(19b) 



Besides this, if the summation index n takes the same 
values in both expansions, the restrictions on the param- 
eters, the form of the recurrence relations and the order 
of the Bessel functions are the same for both sets of solu- 
tions related by (|19a[) . In this event both solutions have 
same characteristic equation and their series coefficients 
are proportional to each other. However, in some cases 
one series breaks off on the left while the other termi- 
nates on the right and so these statements are not true, 
as in the example considered in section III. 

Next we write the first and the fifth set of solutions, 
specify the restrictions on the parameters and the con- 
ditions to use each of the preceding recurrence relations 
(the other sets are in Appendix B). Thus, 

oc 

U[^\z) = E(-l)"^^''^2lU-i (2^/gi) , 

n=0 



U^,^\z) = (z-z,) 



(l-B2)/2 



M.-i (2 V9(^ - ^0)) , (20a) 
[B2^0,-l,-2,---]. 



For U^^ we have [li 



v(l) 



qzo 



{n + 1) 



Bi 



(n+f ) (n+f + i) 



Pi^^ = AB3 - 2qzo + An {n + B2 - I) 
2gzo(f -1) (f + f 



(n+^-l) (n+f ) 



7'- 

In 



qzo{n + B2-2)(n + B2 + ^-l] 



(20b) 



in the recurrence relations for bn '' which are given by 

Eq. (HH) if B2 7^ 1, 2; Eq. (UnD if S2 = 1; 

Eq. dm) if B2 = 2. (20c) 

For J/g"'-' we get the following coefficients {Pn' = Pn') 



qz,in+l)[n + B2 + ^ 



a,, = -- 



Bi 



qzo {n + B2- 2) - 1 - 

h(5) 



in the recurrence relations (|20cp for the &„ . We find the 
formal relation 



^ {~irT[n-{B,/zo)] 
" T[n + B2 + {B,/zo)] " ■ 



(21) 



Actually, this relation holds only if (— Bi/zg) and (B2 + 
Bi/ zq) are not zero or negative integers; on the contrary, 
one solution truncates on the left (a„ = for some n) 
and the other on the right (7^ — 0). There arc similar 
relations for the other solutions. 
From relations (jA.4p . we get 



lim C/f (z) - e^V-J^ z 



2*V9^ Jl/4)-(B2/2) 



lim C/f^z) 



-2i^ (l/4)-(B2/2) 



(22) 



that is, each solution presents one of the possible behav- 
iors given in equation ([5]). On the other side, equation 
([TT|) gives the solution u[^'' as a linear combination of 

Ui'^' and J/}^'' and, consequently, its behavior at z — 00 
must be a linear combination of the previous expressions. 

This implies that the solution J/^^^ is bounded for z = 00 
only if both C/p^ and U^"^' are bounded too. The same 
holds for J/f ^ 

The previous solutions lead to the usual solutions in 
series of Bessel for Mathieu equation. This equation has 
the form 



dv? 



cr2[a- 2fc2cos(2crw)]i(; = 0, g = fc^, (23) 



where a = \ or a ~ i for the Mathieu or modified Math- 
ieu equation, respectively. Then, by introducing 

w{u) — U{z), z = C0S^((TU), 



1 P a 

zo = 1, Bi^ --, S2 = 1, -^3 = y - ^, g = k'^, 



(24) 
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into equation ([5]) we obtain the above Mathieu equation. 
Using these relations we can express the solutions of the 
Mathieu equation in terms of trigonometric (hyperbolic) 
functions. These solutions are given by infinite series 
and the solutions arising from the second subgroup have 
exactly the same characteristic equations as the ones re- 
sulting from the first subgroup. For example, 



w 



(i) 



= V(-l)"6Wz«[2fccos(a^)] 



n=0 



U) 



H = V6«z(^^[2A:isin(a^)], 



n=0 



(25a) 



with the recurrence relations 
1.(1) .^,(1) _ n 



qbY> - a^y' - 0, qb'^' + [4 - a]b\'^ + 2qb\,'' = 0, 



gfeW 1 + [4n2 - a] 6(1) + qb^^U =0, {n> 2), (25b) 



which follow from the second relation given in (|20cp since 
B2 = 1- Both sets have the same coefhcients because 



equation ((2T|) gives fe„ = (—!)"&„. According to the 
next section, the solutions w^"* converge for | cos(tTu)| > 
1 whereas Wg"'"' converge only for |sin((TM)| > 1 if j = 
2, 3, 4. Despite this, in a common domain of convergence 



is a constant multiple of w^?^ , except if fc = |18l.l2 



For the general case, we have found no proof for the 
equivalence between the sets of solutions of the first and 
second subgroups. In fact, equivalence could take place 
only between two sets of infinite-series solutions where 
the summations begin at n = 0. For instance, for the 
first and fifth sets, this requires that neither (— Bi/zq) 
nor B2 + [Bi/zq) are zero or negative integer since this 
implies infinite series in both sets. 

On the other side, some sets of solutions admit the 
limit 0, in which case equation ([5]) becomes a 

double-confiucnt Heun equation with two (irregular) sin- 
gularities located at 2; = and z = 00. In this case the 
second subgroup is irrelevant and only the sets Ui\z) 
and U^^^ (z) admit such limit, as we can see by examining 
the recurrence relations for the series coefficients. 



extending from negative to positive infinity) by introduc- 
ing a characteristic parameter v into the solutions but, 
then, the domains of convergence are changed. 

For n — > 00 the recurrence relations for fel,^'' become 
(we write 6„ instead of bn^) 



qzo 
[4n(n 



1 / Bi 1 
l--[B-2 + — -- 
n \ Zq Z 



0(4 



bn 



+ 1 



B. 



1) 



0(1)" 



qzo 



1 



1 / Bi 1 

n \ Zq 2 



OK. 



bn-1 



whose minimal solution for zq 7^ is 



Un+l 
bn 

bn-1 



qzp 

4n2 
4ri^ 
qzo 



1 



O 



O 



(26a) 



where the notation having the form /„ ^ gn means that 
gn is an asymptotic approximation to as n — cxd, that 
is, fn/gn = 1 when n ^ 00 [2lj. Similarly, for zq — 



bn+l 

K 

bn-1 
bn 



qBi 


H 






4n3 




4n^ 




FO 






1 - 




Wi 









(26b) 



Firstly we regard the expansions in series of Bessel 
functions of the first kind. Using the first relation given 
in equation (jA.7p . we find 



J: 



2n+B-2 + l 



qz 



qz 



1 + 



J2n+B2-1 ['^^/qz) 

Then, keeping only the leading terms, we obtain 



bn J2n+1 



zpq^z 
16n^ ' 



(^0 ^ 0) 



bn+i J2n+B2+i {'2^/qz) Biq^i 

bn J2n+B2-1 (2^/9^) l^Tl^ 



(zo = 0). 



Therefore, the series converges for any finite z, but the 
ratio test is inconclusive when z ^ oo. For this case 
we proceed as in the case of the Mathieu equation [3 ■ 
Thus, by using equation (|A.4p we write (z —^ oo) 



C. Convergence of the solutions 

Now the ratio test is used to get the convergence of 
the first set of solutions. For the other sets the conver- 
gence follows from the transformation rules. The follow- 
ing convergence regions refer to one-sided series; these 
series hold only if the differential equation has an arbi- 
trary parameter that allows to satisfy the characteristic 
equations. If there is no free parameter, it is possible to 
obtain convergent two-sided series solutions (summation 



1 



^2„+B,-i(2Vg^)~(-i)" 



2^ - ^{B2 - 1)7T - ^TT 



By inserting this into u[^\ we find 

lim [/(i'(z) = z(i/4'-(^^/2) ^ 



cos 



I a.Tgy/qz\ < TT. 



2V^-i(i?2-l)7r-i^ 



Y.biP- (27) 



n=0 
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Since the series J2^=o ^ converges, the solution u[^\z) 
also converges at z = oo. However, if the cosine is ex- 
pressed in terms of exponential functions, we obtain a 
linear combination of asymptotic behaviors of the expan- 
sion in Hankcl functions, as observed after equation ([22| . 
Thus, depending on the values of B2 and | arg(y^5z)|, the 

solution C/}^^ may be bounded or unbounded at z = 00. 

Now we consider the expansions in series of Hankel 
functions and Bcsscl functions of the second kind. For 
a = 2n + B2 — 1 and x = relations (|A.7P give 



zol > |2o| if 



Re 



Re 



> 1 : L/^^^ and U^^\ 
< 1 : Ui^^ and ui^\ 



(29b) 

The above conditions have been derived for infinite series 
where the summation begins at n = 0, as in the case of 
the Mathieu equation. A counterexample occurs in the 
analysis of the inverted potential, as explained after the 
paragraph containing the solutions (j49p . 



7U) 



+ 1 [^V^] 



7U) 

qz 



-1 [2v^] 

, 2^2-1 , B2[B2-1] 



2n 



[j-2,3,4]. 



Combining this with the results given in equations (j26ap 
and dMEl, we find 



bn+iZ. 



2n+B2 + l 



1 



z 



1 



B, 



Bi 



^n+1^2»+B. + l (2Vg^) 
bnZ^AB,_, (2VP) 



El 

nz 



O 



1 + 



, 20 = 0. 



Then, the ratio test implies that the series in the solutions 
U[^'' (j 2,3,4) converge for \z\ > if zq = 0. On the 
other hand, for zq 7^ we have {n 00) 



bnZ, 



U) 



qz) 



l^ol 



l + lRc(i?2-2 + :^Uof4 

20 



, (zo ^ 0). (28) 



III. THE INVERTED POTENTIAL 

In this section we get solutions for equation ([3|) of the 
inverted potential by using the preceding solutions in se- 
ries of Bessel functions of the first kind for the Whittaker- 
Ince limit of the CHE. We will find one pair of finite-series 
solutions and two pairs of infinite-series solutions, all of 
them convergent and bounded for any value of the vari- 
able u. The two pairs of infinite series have the same 
characteristic equation, but only the solutions of one pair 
can be expressed as a linear combination of two series of 
Hankel functions converging for all values of the indepen- 
dent variable. 

The substitutions 



coshul ''*'^C/(z), z 



ip{u) 

bring equation ([3]) to the form 



- I 



sinh li, 



dU 
dz 



(30) 



£ b^ / h\ 

+ z- 1) 

4 16 V2 4/ 16^ ' 



U = 0, (31) 



which is Whittaker-Incc limit ([5]) of the CHE with 

20 = 1, ^1 = -^' ^2 = ^-/, 



Thus, if 2o 7^ in general the scries converge for |z| > 
|zo| since the right side of this equation is < 1 but, by 
the Raabe test [H, [23| , they converge absolutely also for 
|z| = |zo| if the numerator of n is < — 1, that is, if Re(i32 + 
Bi/zo) < 1. This possibihty was not noticed before [T3 | 
because the term of order l/n was not considered. 

The convergence for the other sets of solutions arc 

(i) 

obtained by applying the transformations rules to 1/^ . 
For j = 2,3,4, we find that these solutions converge for 
|z| > |zo| in the first subgroup, and for |z — zo| > |zo| 
in the second subgroup. The special cases are given 
by: 



Re LB 



\z\ > I20I if 



i) < 1 : U[^^ and U^" , 



Rc(B2 + ^) >1: U^^^ and U. 



U). 
4 ' 



(29a) 



B-i = 

4 16 



. (32) 

16 ^ ' 



Thus, the Schrodinger equation ^ will be solved by re- 
placing U{z) in equation (j30p by the the expansions in 
scries of Bessel functions of the first kind with the pa- 
rameters ([5^ . Notice that for b = 0, the potential V{u) 
reduces to a hyperbolic Poschl- Teller potential whose so- 
lutions are given by hypergeometric functions [23 | . 

In section HI. A we write down the pair of finite-series 
solutions which allow to get the 'quasi-solvable' part of 
the energy spectrum. We prove that these are degener- 
ate for any finite value of the parameter I and adapt a 
method devised by Bender and Dunne [26| (for a poten- 
tial which leads to a biconfluent Heun equation) to find 
the energy levels and the series coefficients. In section 
III.B we discuss the infinite-series solutions. In all cases. 



using the parameters ([5^ we find that the order of the 
Bessel functions is half-integer and so these functions are 
represented by finite series rather than by infinite ones 
j25l |. Besides this, the recurrence relations are always 
given by equation p2|) since a_[ = 0. 



where the recurrence relations for the coefficients P„ 
are 



Pn+l + finPn + InPn-l = 0, = 0) (36a) 



with 



A. Finite-series solutions 

According to section II. A, if 7n=Ar+i = in the re- 
currence relations (|12p. the series terminates at n = 
and we obtain a finite-series solution. In this case the 
recurrence relations can be written in the form 



f ao • • • 

71 Pi cei 
72 P2 0L2 



\ ••• 







\ 



JN-I Pn-1 OLN-1 hN~l 

7jv 

= 0. (33) 



This system has nontrivial solutions only if the determi- 
nant of the above tridiagonal matrix vanishes. Besides 
this, if (as in the present problem) the elements of this 
matrix are real and 



a,7.:+i>0, 0<i<A^-l, 



(34) 



then the A^ -I- 1 roots of the determinant are real and 
distinct 

Only UY' and U!^" give finite-series wavefunctions. In 
order to prove the degeneracy and apply the procedure 
of Bender and Dunne, we redefine the series coefficients 
as 



6(1) = 



T[2n + B2] Pn 



{-qzoY n! r[n- (Bi/zo)]' 



(2) ^ r[2n + 2 + E2 + (2Bi/zo)] Qn 
{-qzoY n\T[n + 2+{Bi/za)Y 

Then, inserting the solution u\^^ and into equation 
([50)1 and using the parameters given in we get the 
following pair {ijjf,tpi) of even and odd solutions 



1 /A 
V'?(«)- [tanhi.]'-^5^(-l)" - 



r[2n^l + (3/2)] P„ 
n\r[n + (1/2)] 



-^2ri-i+(i/2) ( -sinhu ) , 



(3n = -Ei-n[n-l + - 



b2(2/ + 1)(2Z- 1) 



32(4?i- 1 - 2l){4n + 3-2l)' 



n{2n-l){2n~2l-l){n-l) 
64 (4n + 1 - 20(471 - 3 - 20(4n - 1 - 2iy ' ^ ' 



E, := 



£ b^ 



I 1 



4 32 V2 4, 
The recurrence relations for Qn are 

Qn+l + PnQn + InQn-l = 0, {Q -i = 0) (37a) 

with 

r 



I3n = -Ei-{n-l + l)\^ + 

h^(2l + l){2l - I) 
32(4n+ 1-20(471 + 5 -20' 



In 



b^ 77(277-f l)(2n-2/-hl)(n-/) 
64 (477 - 1 - 20(477 + 3 - 2Z)(4n + 1 - 20^ 



(37b) 

From the relation J\{—x) = {—l)^Jx{x), it turns out 
that these solutions are even and odd, that is, ipf{—u) = 
^P1{u) and 7A?(-m) = -i^iiu). 

The series are finite because the coefficients 7„ and 7„ 
of Pn-i and Qn-i vanish for tt, = I and, consequently, 
the series terminate at 77 = Z — 1 as stated above. In 
virtue of equation p4[) . each cigcnfunction corresponds to 
/ distinct and real eigenvalues. In addition, from J\{x) — 
(x/2)^/r(A + 1) when x ^ 0, wc find 



lim ?/'i(?7) ~ [coshw] 



finite, 



and from the first of equations (|A.4[) we get 



u — >±oo 



[tanh-i 



V sinh u 



■ COS 



— sinh 7i 
2 



r[277 -? + 3/2]P„ 

77! r [77 + 1/2] 



ipl{iL) = [tanhit]' ^ 



i-i 



ri=0 



2n 



r[277 -Z + (5/2)]( 



77! r [77 + (3/2)] '-J^^-^Hm)[^^^^'-^]^ii^) 



Thence, tpl{u) is bounded also at the singular points of 
the equation (z = and z = —00). The same holds for 
the solutions ipKu). In fact, expansions in series of the 
other Bessel functions are also bounded for all values of 
u but do not present definite parity. 
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The degeneracy of the previous solutions is estabhshed 
by arranging the recurrence relations in the matrix form 
([55)1 . For Pn we write AP = 0, where 

P=(Po,Pi,... ,P,_i)* 

(t means 'transpose') and 



/ Po 1 
71 /3i 1 

72 ^2 



V 



\ 



11-2 Pl-2 1 

Pi-i J 



(38) 



For Qn wc write BQ = 0, where 



and 



//3o 1 
71 /3i 1 

72 /32 



V 



\ 



7;_2 A-2 ^1 

7;_1 / 



This matrix can be rewritten as 

/ Pi-i 1 
7/_i /3/_2 1 

7i-2 Pi-z ■ 



\ 



72 Pi 1 
71/^0/ 



(39) 



V 

due to the identities 

/3o = /3i = A-2, /32 = A-3, • • • , A-1 = /3o, 

7i=7/-ii 72=7;-2, 73=7i-3, 7/-i=7i- 

Therefore, the two matrices are constituted by the same 
elements. To prove the degeneracy, it is sufficient to 
show that these matrices possess the same roots, that 
is, detA = detB. This is obvious for ^ = f and 1 = 2. 
For I > 3, we use the ^-by-Z antidiagonal matrix § hav- 
ing I's on the antidiagonal as the only nonzero elements, 
that is, 



1 



detS : 



1 



Then we find the similarity relation 

A = §-iB*§ = §B*§, 



(40) 



(41) 



where B* is the transpose of B. Thus, from the prop- 
erties of the determinants it follows that detA = detB 



and, therefore, the finite-series solutions are degenerate 
for any finite /. 

On the other hand, the eigenvalues may be computed 
by equating to zero the determinants of the preceding 
matrices. However, the procedure of Bender and Dunne 
gives as well the coefficients P„ and Qn as polynomials 
of degree n in the the parameter Ei . The procedure is 
implemented by taking Pq = Qq = 1 as initial condi- 
tions and by using the recurrence relations to generate 
the other coefficients. For a fixed I, the eigenvalues are 
obtained by requiring that P; = or Qi = 0, since the 
series terminate at n = i — 1. Thus, equations (|36ap and 
(l36b)l yield 

P„+l = [El + knl]Pn - 7n^'n-l, P-l = 0, Pq = 1 (42a) 



with 



kni = n\ n-l + 



wherefrom we find 



b^(2Z + l)(2Z-l) 
2) 32(471 - 1 - 20 (4n + 3- 2/)' 

(42b) 



Po = 1, Pi = Ei + ko 



P2 = Ef + [koi + ku]Ei + koiku - 71, 



P3 = Ef + [koi + kii + k2i]Ef + [koikii + koik2i 
+kuk2i - 71 - 72]^'; + k2i [koiku - 71] - fco/72, (43) 



and so on. Similarly, we write the relations (|37ap and 
(|37bp for Qn as 

Qn+l = [El+knl]Qn-lnQn-l, Q-1 = 0, Qo = 1, (44a) 

where 



[n-l + 1] 



"+2 



(21 + 1) (21-1) 



32(4n+l-20(4n + 5-2/)' 
(44b) 

Then, the expressions for Q„ are obtained by replacing 
kni and 7„ by k^i and 7„ in the expressions for P„. 

As an example we find the energies and the respective 
eigenfunctions for I = 1 and 1 = 2. For / = 1 the energy 
that follows from the condition Pi = Qi = is 



Ei-_ 



32 



= 



£ = - (b^ - 1) , (45a) 



corresponding to the degenerate pair of eigenfunctions 



— sinh u 
2 



(u) = V tanh u J_ 



ipi(u) = V tanh u Ji ( — sinhw 



(45b) 
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For / = 2 the condition P2 = Q2 = leads to 



^/=2 - 4 



£± = i [b^ - 5] ± ym?. 



(46a) 



Since, Pi = E2 + {ih^ /Z2) and Qi £;2-(l/2)-(5bV32) 
these energies are associated with the eigenstates 



= (tanhii)' 



J_3 



■ sinh -i 



4- f — + 1 ± Vb^ + 1^ Ji ( - sinh u 
bM 2 / ^ V 2 



= (tanhu)^ 



■ sinh u. 



1 I J3 I — sinh u 



, (46b) 



up to normahzation factors. 

These energies for / = 1,2 arc the same found by 
Cho and Ho [5|. However, the cigcnfunctions differ from 
theirs, as we can see by writing the Bessel functions in 
terms of elementary functions via the formulas written 
in Appendix A. Indeed, their solutions may be obtained 
from another group of expansions in series of Bessel func- 
tions given in [l^ . 

Notice that the coefficients of the preceding finite series 
factorizc in the same manner as the coefficients of the 
problem considered by Bender and Dunne (26j . that is, 
for a fixed / one has 



Pi+,{Ei)^p,{Ei)Pi{Ei), 
Qi+,{Ei)=q,{Ei)Qi{Ei), (i>0) 



(47) 



where Pi and qi are polynomials of degree i in £';. For 
example, by taking n = / (7; = 0) and n = Z + 1 in 
equation (|42ap . we find 

Pz+i = {Ei+ki^i)Pi ^piPu 

Pi+2 = [{El + ki+u){Ei + ki^) - 7i+i]Pi = P2Ph 

and; by induction, wc obtain the previous expression for 
Pi+i- 



B. Infinite-series solutions 

The solutions U^'' and U^'' lead respectively to odd 
and even infinite-series wavefunctions which are bounded 
for any value of u, namely, 

00 

V'2°(u)= [tanhw]-'-^^(-l)"6(f) ^ 



Tl = 



J. 



2n+/+(3/2) 



■ sinh u 



n=0 



J2„+;+(i/2) ( ^ sinhii ) . (48a) 



In the recurrence relations (fT2|) for hn^ the coefficients 
are 



a, 



(3) 



b^(?i + l)(2?i + 3) 
' 2{'in + 2l + b) (4n + 2/ + 7)' 



+ 



1: 



(3) 



\? (2/ + l)(2/-l) 
;(4?i + 2; + l)(4n + 2/ + 5)' 



b2(2n + 2Z + l) [n + l) 
'2{An + 2l- I) {An + 21 + 1)' 



(48b) 



and for Ihif' the coefficients are 

b2(n-|- l)(2n-t- 1) 



a. 



(4) 



2{An + 2l + i){An + 2l + hy 



?(4) 



+ I 



+ A{n + l)[n + 



1 

8 ■ 2 
b^ {21 + 1) {21 - 1) 
; {An + 21-1) {An + 21 + ?,)' 



1: 



(4) 



b2(2n + 2/-l) {n + l) 
2{An + 2l-i) {An + 21- I)' 



(48c) 



According to section 2.3, the corresponding expansions 

fi 21 

given by series of the Bessel functions and Ha ' con- 
verge in the domain \z\ = sinh^ m > 1 which does not 
include all the values of u. 



On the other side, from the solutions Uj^^ and U^''' we 
find another pair of bounded eigenfunctions given by 

rs{u) = tanh.f:Il|^^(-l)"6(f) 



r(i) 



ri=0 



X J. 



2n+; + (3/2) 



■ cosh u 



X J. 



2n+i+(l/2) 



■ cosh u 



(49) 



in which the series coefficients are proportional to the co- 
efficients of the previous pair and, consequently, we have 
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the same characteristic equations. Since Re(i?i/zo) < 1 
in equation (j29bp . now the expansions in terms of Yq, and 

i/i^'^' converge in the domain |1 — = cosh^ m > 1 which 
covers the entire range of u. 

By using the solutions ui^^^ and J/g"''' we would find so- 
lutions equivalent the preceding ones. Actually, wc would 
find an''^^ cc {n ^ I + 1) in the recurrence relations for 
bn'^ and bn"^ what means that the series begin at n = /. 
Thence, by setting n = m + ^ we may conclude that such 
solutions are proportional to the above ones. 

At last, notice that the previous considerations take 
into account only the analytical properties of the wave- 
functions. The full solution of the problem requires the 
computation of the characteristic equation resulting from 
the three-term recurrence relations, which is represented 
by an infinite continued fraction having the form given in 
equation (|12p or by the determinant of an infinite tridi- 
agonal matrix. 



IV. THE ELECTRON IN THE FIELD OF A 
POINT DIPOLE 

For an electron with mass M, charge e and energy E in 
the field of a point electric dipolc, the time-independent 
Schrodinger equation is 0] 



2M 



D ■ f 



which is the Whittaker-Ince limit ([5]) of the CHE, with 
the following set of parameters 

zo = 1, Bi = -m-l, B2 = 2m + 2, 



-2/3. (53) 



- S U. = 0, 



Therefore, the solutions 8(6') can be constructed by in- 
troducing solutions of equation ([5]) into equation (|5ip . 
Then, the admissible values for C are determined from 
the characteristic equations which follow from the recur- 
rence relations for the series coefficients. 

For f3 — the angular equation (|50ap has solutions 
regular in the interval 0<6'<7rifC = £{£ +1), where 
^ is a non-negative integer such that i > m. These solu- 
tions are given by the associated Legendre polynomials 
Pp{cos9). For this trivial case, a closed form for energy 
spectrum can be obtained from boundary conditions on 
the radial part of the wavefunction. However, ii (3 ^ 
there no analytic formula for £ since the constant C must 
be determined from a transcendental equation (charac- 
teristic equation). 

In section IV. A we write the expansions in series of 
Bcssel functions for the angular equation. Wc find two 
periodic expansions having period 27r and the same char- 
acteristic equation. However, in section IV. B we find 
only one expansion in series of associated Legendre poly- 
nomial; this has the same characteristic equation as the 
solutions in terms of Bessel functions. 



where D is the dipole moment. This equation is separable 
in spherical coordinates {r,0,ip). Choosing the z axis 
along the dipole moment and performing the separation 



1 



'^ir,e,ip) ^ -R{r} 9(6') e 
r 

{0 < < tt) one gets 



0,1,2,... 



1 


d 










sin 6 


Me 





sin^ e 



R = 0, 



9 = 0, 
(50a) 

(50b) 



where C is the separation constant, /3 = 2MeD/h^ and 
£ = 2ME/1t?'. The energies are determined from the 
solutions of the radial equation, but firstly it is necessary 
to determine the parameter C from the solutions of the 
angular equation. The substitutions 



9(6') = (sin 61)"' U{z), z = sin^ (6'/2) , 
jive the equation 

.<PU , , , ,dU 



(51) 



dz 



[m{m+l)-P-C-2f3{z-l)]U = Q, (52) 



A. Expansions in series of Bessel functions 

For P the solutions 9(0) can be obtained by in- 
serting into ([51]) the solutions in series of Bessel functions 
of the first. 

The other expansions are unsuitable in virtue of their 
domain of convergence. 

Thus, using the parameters ([53]) we find out that only 



U[^\ U^'^\ U^'^'^ and U^'^'^ are vahd. Moreover, both u[^^ 
and [/j^-* afford the same solution, denoted by 9i. Anal- 
ogously, both U^^'' and Uq^'' yield another solution, de- 
noted by 92. These solutions are 



e 



-1 / ^\ m oo 

em = (sin-l cot- ^(-1)"6W 

/ \ / n=0 



X J2«+2m+i ( 1^8/3 sin - ) , (54) 



em = cos^ tan^ Y.biPx 

J2n+2m+l (V&P ^OS M , (55) 
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with the recurrence relations for bn^ given by (a_i = 0) 



functions in terms of Jacobi's polynomials Pn"'^^ through 
equation (jA.Sp . we find 



P{n+l) ,(1) 



(2n + 2m + 3) 
m(m + 1) -C 



-'n+l 



b^ 



ri(n + 2m + l) + 

(3{n + 2m) ^(i) 
(2n + 2TO- 1) ""^ 



0. (56) 



These solutions have the same coefHcients bn^ because 
equation ([2T|l implies that bn^ = (— l)"'6n^ Thus, there 
is only one characteristic equation to determine the val- 
ues of the parameter C. 

The solutions Oi and O2 are connected by the substi- 
tutions 6 ^ 6 + t: and (3-^-/3 which leave invariant 
the angular equation (notice that this is equivalent to 
the change of &„ by (— l)"6,i in the recurrence relations). 
Both are convergent in the interval < 6 < it and are 
regular at the singular points 9 = and 9 = tt. We 
have found no criterion to discard one of these solutions, 
neither have found a proof that they are equivalent (the 
problem mentioned in section 2.2). However, in the fol- 
lowing we will find only one solution in series of Jacobi 
polynomials; this has the same characteristic equation as 
the above solutions. 



B. Expansions in series of Jacobi polynomials 

The expansions in series of Bessel functions given in 
section II appear associated with a group of expansion 
in series of Gauss hypergcomctric functions F{a,b; c; x), 
which arc obtained by applying the transformations rules 
to the solution 



Ui(z) 



00 

n=0 



-n,n + B2 - l;B2 



Bi 



(57) 

[B2 ^ 0, -1, -2, • • •] where is the same as in (j20ap 
and satisfies the recurrence relations (|20cp . This series 
converges for finite values of z. The restrictions on the 
values of B2 assure independence of the hypergeometric 
functions. However, we must demand as well that c ^ 
0, — 1,— 2,--- because in general F{a,b;c;x) is not de- 
fined if c is a negative integer or zero. For this reason, 
this group is less general than the group formed by series 
of Bessel functions. Now, by taking 115(2;) = 5^Ui(z), 
we find the solution 



00 

U5(z) = ^ {-n,n + B2- 1; - — ; - ) , (58) 

n=0 



Bi z 

zq' zo^ 

[B2 ^0,-l,-2,...] 



(5) 

where the coefficients 6„ are formally connected with 
blP by equation pT|) . If B2 + Bi/zq and (—Bi/zq) are 
not zero or negative integers, both solutions are valid and 
are given by infinite series. Then, setting z = zq cos^(crit), 
using the relation (pij) and rewriting the hypergeometric 



Ui(z)=Wi(ii) = T{a + 1)Y^ 



^l blP 



r{n + a + 1) 



n=0 



X P^"''3)[c0s(2crM)], 



U5(2)=W5(w) 



r(/3 + i)^ 



(-1)" nl b'^^ 
r(,i + a + l) 

cos(2(Tw)] , 



where a ~ B2 — 1 + Bi/zq and /? = — 1 — Si/zq (this /3 
should not be confused with the parameter of the angular 
equation). Then, relation (jA.lOp implies that Ui is a 
multiple of U5. This conclusion does not hold if only 
one solution is valid. The same can be said of the other 
pairs of solutions, (Ui, Ui-1-4). In this manner, the linear 
dependence of solutions having the same characteristic 
equation is almost trivial in the present case. 

The solution Ui (equivalent to IU5) is the only one 
applicable to the angular equation of the point dipole. 
Thus, by inserting this into l|5ip and using the parame- 
ters dSni), we find 

00 

m = (sine)" ^6(1) 

n=0 

X F[-n, n + 2m. + 1;to+1;cos^(6'/2)], 

where the coefficients bn^ again satisfy the recurrence 
relations ([55]) . In terms of Jacobi polynomials (jA.8p or 
associated Legendre polynomials (|A.11[) . we find that 



(sine)-V)zi^_6(i)p(™.n)(, 



{n + my. 



COSf 



^ — ^ I TJ A- / m 1 1 



^ (n + 2m) 



(59) 



If m = 0, this reduces to the solutions in series of ordinary 
Legendre polynomials Pn ~ Pn given by Levy-Leblond 
0. If m ^ 0, these are the Alhaidari solutions in terms of 
Jacobi polynomials [7| up to a redefinition of the series 
coefficients. On the other side, by taking into account 
the solutions U2, U3 and U4, written in [iJIj we find 
that these are not valid by one of the following reasons: 
the last parameter of the hypergeometric functions is a 
negative integer or 8(6') is not regular at 6* = or = tt. 

Finally, the expansions in hypergeometric functions are 
inapplicable to the Schrodinger equation with the in- 
verted potential In effect, infinite-series solutions 
are inappropriate because they converge only for finite 
values of the arguments of the hypergeometric functions. 
On the other side, only Ui(z) and 112(2;) could afford 
finite-series solutions but in these cases the hypergeomet- 
ric functions are not defined since the last parameter is 
zero or negative integer, that is, c = B2 + Bi/zq = 1 — 1. 
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V. CONCLUSION 

We have dealt with solutions for the Whittaker-Ince 
limit (O of the confluent Heun equation (CHE) and pos- 
sible applications of these solutions. Specifically, in sec- 
tion II we have considered expansions in series of Bessel 
functions of the first and second kind, in addition to so- 
lutions in series of Hankel functions given in a previous 
paper jT^]. In sections III and IV we have established 
solutions in series of Bessel functions for an inverted po- 
tential and for an angular equation for a point electric 
dipole, respectively. 

Under certain conditions the expansions in terms of 
Bessel functions of the first kind are convergent and 
bounded for any value of z, including the point z = oo. 
These are suitable to solve angular equations since in 
this case z assumes finite values around z = 0. The 
other expansions in general converge for \z\ > \zq\ or 
\z — zqI > l^ol but, under special conditions, may con- 
verge also at \z\ ~ |zo| or |z — zo| = |zo|- For this reason, 
these and the solutions in series of Bessel functions of the 
first kind are applicable to the inverted potential. 

We have noticed the possible coexistence of different 
sets of expansions in series of Bessel functions having the 
same characteristic equation. For the Mathieu equation 
it is known that these solutions are linearly dependent in 
a common region of convergence [3 ■ The same holds 
for expansions in series of Jacobi polynomials, as we have 
shown in section IV. B. If this dependence is valid for the 
general case, we can avoid the duplicity of solutions in 
infinite series of Bessel functions. Nevertheless, we have 
found no proof for such conjecture. 

In section III wc have obtained eigenstates given by 
finite and by infinite series of Bessel functions for the 
one-dimensional Schrodinger equation with the inverted 
quasi-exactly solvable potential. For any value of the 
parameter I which characterizes the quasi-exact solvabil- 
ity, we have proved that the degeneracy of even and odd 
states given by finite series follows from the three-term re- 
currence relations for the series coefficients. These quasi- 
polynomial solutions permit to determine only a part of 
the energy spectrum by using, for instance, the procedure 
of Bender and Dunne presented in section III. A. 

We have found that the infinite-scries solutions are also 
convergent and bounded for any value of the variable u. 
Odd and even solutions are given by the expansions in 
series of Bessel functions of the first kind, while other 
kinds of Bessel functions give bounded solutions without 
definite parity. In principle, the characteristic equations 
associated with these solutions may afford the remaining 
part of the energy spectrum. 

In section IV we have found solutions given either by 
series of Bessel functions or by scries of Jacobi polyno- 
mials for the 6'-depcndencc of the scattering of electrons 
by the field of the point electric dipole. The expansion 
in series of Jacobi polynomials is equivalent to the one 
found by Alhaidari [7[ and includes, as a particular case, 
the solution in series of Legendre polynomials proposed 



by Levi-Leblond for the case m = [Gj]. We have found 
two expansions in series of Bessel functions correspond- 
ing to the same characteristic equation. We have also 
seen that the expansions in series of Jacobi polynomials 
arc inappropriate for the inverted potential. 

We observe that, by using transformations rules, we 
can find different sets of solutions with the same charac- 
teristic equation also for expansions in series of confiuent 
hypergeometric functions for the CHE, given in [l^l - 
these are the ones which lead to the expansions in series 
of Bessel functions discussed in section II. Then, the is- 
sue concerning the linear dependence or independence of 
solutions arises in this case as well. 

Finally we remark that we have dealt only with one- 
sided series solutions which are valid if the differential 
equation has an arbitrary parameter. If there is no free 
parameter, we have to deal with two-sided infinite series 
solutions (— oo < n < oo). For the CHE there are, for 
instance, expansions in double-infinite series of hyper- 
geometric and confiuent hypergeometric functions 0, [131 
which admit of the Whittaker-Ince limit [ij] . These solu- 
tions need to be extended to incorporate all the Meixner 
expansions in scries of Legendre and Bessel functions for 
the ordinary spheroidal wave equation [28l . [29j . that is, 
for the CHE © with = 0. 



APPENDIX A: SOME MATHEMATICAL 
FORMULAS 

Firstly we give some properties of the Bessel functions 
Ja{x), Ya{x), Ha\x) and Ha\x). After this we write 
some formulas concerning the hypergeometric functions. 
The power-series representation for Ja{x) is 



(-1)' 



k=0 



k\V{a + k + l)\2 



2k 



(A.l) 



This function is connected with the Hankel functions by 
the relation 



Ux) = \[H^^\x)+H^^\x)]. 



(A.2) 



Similarly, the Bessel functions of the second kind can be 
expressed as 

Y^{x) = Uh^^\^)-H^I\^)\. (A.3) 
Zi 

On the other hand, for a fixed a the behaviors of the 
Bessel functions when |a;| -^■ oo arc given by psi ] 



Ja{x) ~ \ —— COS ix — iftTT — ^TT 1 , | arg x\ < tt; 

\ TTX \ 2 4 



-TT < arg X < 27r; 
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H^aH^) ^ \/— e-*(='-5«'^-3'^), -27r < arg x < tt. 
V nx 

. (^-^^ 

The behavior for Ya(x) is obtained by changing the 
cosines by sines in the expression for Jq(x). The Bessel 
functions are given by finite series of elementary func- 
tions if their order a is half an odd integer; then the 
restrictions on argx are unnecessary in (|A.4p p^ . For 
instance, if m = 0, 1, 2, • • • , the Bessel functions of the 
first kind can be expressed as 



J_ 



Ax) 



d \ cos X 
xdx I X ' 



which give, in particular 
'2 



J- 1 (x) 



COSX, Jl (x) 



sma;, 



2 / . cosa;^ 

sin a; 
nx \ X 



sm X 



(A.6) 



where n is a non-negative integer. The Jacobi polynomi- 
als can be expressed as 



(-1)' 
2"n! 



(l-a;)-"(l+x)- 



., ^- [(l-x)"+"(l + x)'5+"] , (A.9) 

whereby we find the relation 

Pi^^^\-x) ^ (-1)"P("'«(X), (A.IO) 

used in section 4.2. 

On the other side, the relation between hypergeomet- 
ric functions and the associated Legendre functions P^, 
when fc is a positive integer or zero, is given by 



F{k ~ iy,iy + k + l;k + l;y) 



{-l)''k\r{iy + l-k) 



r(i/ + fc + 1) 

X (y-y2)-4p^fc(i_2y). (A.ll) 



Thence, by setting a — p = k in (jA.8[l and ly — n + k in 
(|A.11[) we obtain 



(-2)'=(n + fc)! 



{n + 2ky. ^ 
where we have put ^ = 1 — 2y 



(A.12) 



Since, for x fixed and a — > cx), 
1 



Ja{x) 



r(a + l) \2) ' 



Y^{x) ^ ^tHi'\x) ^ ^H<i\x) ^ --V{a) (-X , 

TT \X J 



we find the relations 
Ja+2{x) 



Ja{x) A{a + l){a + 2)' 



APPENDIX B: THE OTHER SOLUTIONS FOR 
EQUATION (HI) 



From the solutions (|20a[) . the others are obtained by 
using the transformations ^ and as indicated in the 
sequences pT|) and p^ . In this manner, we get 

oo 



n=0 



^ ^2n+l + B2 + (2Bi/zo) C^V^) ' 



Zi'\x) 

which have been used in section II. C. 

In section IV, the relation between hypergeomctric 

functions and the Jacobi polynomials P,!"''^'' was ob- 
tained from [231 



F{-n, n + l + a + l3;l + a;y) = 
nlT{a + l) 
r(n + a + l)^" 



^ -N 1 1 1 -Bi B2 , 



n=0 



(-l)"&(f) 



X Z. 



U) 

2n+l+B2+(2Bi/zo) 



(2^q{z - zo)) , (B.l) 



B2 Bi , 3 

f + i^-l'-2-^'' 



where 



0' 
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/3„ ==4B3-2gzo + 4 



n+ 1 + — 

^0 



n + B2 + — 



2,.o (f - 1) (f + f 



^(3) 
In 



(n+i-f) (n+l-f) 



(B.6) 



(3) 

in the recurrence relations for 6n which are given by 
Eq. (dm if ^2 ^ 2, 3; Eq. ([T31) if B2 = 3; 



qzo 



za 



(B.2) 

(2) 

in the recurrence relations for 6,1 , which are given by 
the equations 

Eq. ©iff + ^^0,-i; 

Z Zq Z 



©iff + ^ = 4 Cliff + ^ = 0. (B.3) 

Z Zq Z Z Zq 

These relations also hold for the coefficients h^n' with 
Hf^ = and 



v(6) 



(^+1)( 



n + B2 + #^ 



Eq. dUl) if S2 = 2. 
We obtain the coefficients /3„ = /?« and 



a, 



(7) 



qzo (n + 1) n + 2-B2 



(n + 2-f)(n+|-f) ' 
gzo (n + 2-B2)(n+l + |L) 



in the recurrence relations (jB.7p for h^r? ■ Finally. 



(B.7) 



(B.8) 



!l , B2 1 



?i=0 



gzo ( ?^ + S2 + 2^ ) ( n + 1 + fJ- 

'"-i + f + f)("+f + f)' 

in the recurrence relations (jB.Sp for the coefficients bi? . 
For the third and seventh sets we find 



1--B2 



^0 z =0 



B2 1 



ri=0 



n=0 



X Z, 



0) 

2n+3-B2 



(2^/9(2-20)) , (B.5) 
[^2^ 4,5,6, •••] 



with the coefficients 

gzo (n + 1) (71 + 2 + ^) 



v(3) 



(n + 2-f) (n+l-f)' 

4^3 - 29Z0 + 4(n + l)(n + 2 - B2) 
2g^o (f - 1) (f + f ) 



T1=0 



^2n+l-B2-(2Bi/2o; 



(2^/9(2-^0)) , (B.9) 



52 Si 3 5 ^ 

2 + zo ^ '2' '2' 



with 



v(4) 



qza 



/^i^) ^ 4B3 - 2gzo + 4 



n - B2 + 1 - — 



2gzo (f - 1) + f 



v(4) 



n + 1 - 



i?2-2^ 



1 _ ^ _ ^ 

2 2 Zo 



-B2 ^1 

2 zo 



(B.IO) 



in the following recurrence relations 



Eq. ©iff + ^^0,i; ©iff + ^ = i; 

/ Zo ^ I Zq Z 
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The coefficients 6„ also satisfy tlie recurrence relations 
(iRTTj) with /3i®' = and 



gzo (n + 1) (?^ + 2-B2 - f- 



1 Bi^ \ / „ I 3 _B2 Bj^ 

^ 2 20 M ^ 2 2 zo 



2 2 zo / V 2 20 
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